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ABSTRACT 



We study the relation between the directed polymer and the directed 
percolation models, for the case of a disordered energy landscape where the 
energies are taken from bimodal distribution. We find that at the critical 
concentration of the directed percolation, the directed polymer undergoes a 
transition from the directed polymer universality class to the directed perco- 
lation universality class. We also find that directed percolation clusters affect 
the characterisrics of the directed polymer below the critical concentration. 
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Several models have been developed for directed path in disordered media. 
Two well known models are directed polymer and directed percolation [1,2]. 
While the directed polymer is based on global optimization, the directed 
percolation can be described as a local process. The directed polymer is 
considered as a model for physical processes such as tearing or cracks [3], 
while the directed percolation has been used to model for example invasion 
of low viscosity fluid in high viscosity [4], as well as the interface of liquids in 
materials such as paper [5]. The relation between the two models is not yet 
clear and even controversial, as discussed bellow. 

The two models of directed polymer and directed percolation can be de- 
scribed in a similar way as follows: In a square lattice which was cut along 
its diagonal and oriented as a triangle whose apex is up and the diagonal is 
its base, we assign to each bond a random number. Of all the paths leading 
from the apex (the origin) to the base we refer only to those whose direction 
is always down to the base. For each one of these paths we calculate the sum 
of the random numbers along it. 

In the directed polymer model, the sum of the random numbers along 
the path defines its value. We focus our interest on the path of minimal 
value (which we call "the optimal path"), and define the roughness exponent 
v by D ~ t", where D is the mean distance of the endpoint of the optimal 
path from the center of the base, t is the size of the triangle, and v is the 
roughness exponent. In this model the random numbers are usually taken 
from a continuous distribution, and thus the optimal path is uniquely defined. 
Huse and Henley [1] who introduced this model found that v = 2/3 and the 
value 2/3 is considered to be exact. Another exponent which characterizes 
this model is u, defined by a E ~ t w where E is the the value of the optimal 
path, and oe is its standard deviation. These two exponents are related by 
the scaling relation u = 2u — 1 [1]. 

In the directed percolation model the random numbers are taken from 
a bimodal (0,1) distribution, and we can define a percolation cluster as the 
collection of lattice sites which are connected by zero sum paths to the origin. 
In this model we can also define a rougheness exponent by W ~ L v where 
W is the mean width and L is the mean length of a percolation cluster. 
W and L depend on the probability to get (denoted by p) by the relations 
W ~ (p c — p)~ u± , L ~ (p c — p)~ Ul] , where p c is the critical probability ~ 0.6447, 
v± is the transverse exponent ~ 1.097, and u» is the longitudinal exponent 
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~ 1.733 [4]. Since v = v±/v\\ — 0.633, there are two independent exponents 
in this model, compared to only one in the directed polymer model. 

Husc and Henley [1] refer to the possibility that in the directed polymer 
model the random numbers are taken from a bimodal probability distribu- 
tion, and stated that the roughness exponent has the same value as in the 
continuous case, i.e. v = 2/3. On the other hand, more recently, Lebedev 
and Zhang [6] claimed that in the bimodal case, for all p < p c the exponent 
is actually the same as in the directed percolation case, i.e. v ~ 0.633. 

In this Letter we address this controversial issue by studying the directed 
polymer model with bimodal (0,1) distribution. The numerical simulations 
lead to the following main conclusions: 

1. For p < p c and for p > p c , for large t the value of the roughness 
exponent is ~ 2/3, the directed polymer value. 

2. For p = p c , the value of the roughness exponent is ~ 0.633, the directed 
percolation value. 

3. For p < p c and t L, the sum of the random numbers along the 
optimal path is inversely proportional to the mean directed percolation 
cluster length, i.e. E ~ t/L. 

It should be emphasized that a study of the directed polymer model yields 
dependence on both the critical probability and mean percolation cluster 
length, L, of directed percolation. So the general conclusion is that the two 
models are closely related. 

In this study the random numbers are taken from a bimodal (0,1) dis- 
tribution, and thus there is usually more than one optimal path and more 
than one optimal endpoint. If there exists a real percolation cluster spanning 
from the origin to the base, all the optimal endpoints (and no other point of 
the base) belong to this percolation cluster, but otherwise, the optimal paths 
and endpoints are not necessarily connected in any special manner. The set 
of optimal endpoints can be characterized by two variables: 

1. Dq - the distance of the center of the optimal endpoints set from the 
center of the base. The center of the set is computed by (Xi + X r )/2 
where X[ and X r are the leftmost and rightmost points of the set. 
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2. W - the width of the optimal endpoints set, which is computed by 
W = (X r -Xi). 

It is expected that for large t, D ~ t VD , and W ~ t Uw , where v D and 
vw are in the range [0,1]. In order to test this hypothesis and estimate vd 
and vw, we define the local roughness exponent of a variable V at the point 
t, as the slope of the curve relating the logarithm of V to the logarithm of t. 
The numerical simulations provide estimates for the local exponents of Dq 
and Wq by computing (\ogV(2t) — log V(t/2))/ log 4, where V is either D 
or Wq. These estimates are presented in Figures 1 and 2 for 5 probabilities: 
p = 0.25, p = 0.5, p = 0.62, p = 0.64 and p = 0.6447(~ p c ). In Figure 
1 it is seen that the slope of the p ~ p c curve is stable at a value ~ 0.63, 
while as p is further below p c , the slope seems to crossover from ~ 0.63 
to a value close to 2/3. In Figure 2 it is seen that the slope of the p ~ 
p c curve converges to a value ~ 0.63, while for p < p c the slope declines 
and approaches (for the probabilities further from p c ) a value close to 1/3. 
These crossovers from directed percolation exponents to directed polymer 
exponents can be explained due to the relation between the triangle size t 
and the longitudinal correlation length £y [2]. For t < £y we expect properties 
of directed percolation, while for t > £y we expect properties of directed 
polymer. 

The results presented in Figure 2 deserve some explanation, as one might 
expect that Wq should be proportional to W - the width of a typical directed 
percolation cluster (in the relevant probability). If this was the case, we would 
get in Figure 2 an initial slope of 0.63, and than the slope would decline to 
zero. Evidently, the picture is different, as for each probability p < p c the 
slope will eventually converge to a value close to 1/3. The reason for this 
lies in the ultrametric tree structure of the directed polymer problem. It was 
shown [7] that in the continuous distribution case, the region of endpoints 
of paths whose value difference from the value of the optimal path is smaller 
than a (small) constant, increases as t 1 ^ 3 . In our case of discrete values, the 
same should hold for zero difference between the values of the optimal paths. 

The main conclusion from Figures 1,2 is that for p < p c and sufficiently 
large t, the optimal endpoints set of width ~ t 1 ^ 3 is located around a point 
whose distance from the center is ~ t 2 ^ 3 . As for t 3> 1, t 1 / 3 is negligible 
compared to t 2//3 , the situation is similar to the continuous distribution case, 
where instead of one endpoint there is a " cloud" of endpoints whose distance 
from the center is ~ 
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For p > p c there is a finite probability for points of the base to belong to 
a real directed percolation cluster, so that Wq ~ t. On the other hand, we 
find numerically that Dq ~ t 1 ^ 2 . As Wo Dq, it is certain that members of 
the optimal endpoints set are found in both sides of the center, and another 
definition of the distance is needed. A definition that takes into account 
the fact that the directed polymer is equally likely to choose any one of the 
optimal paths, is D r = J2i n i\ x i\/ Hi n ii where rii is the number of optimal 
paths whose endpoint is X{. The local exponents of D r are presented in 
Figure 3 for 3 probabilities: p = 0.5, p = 0.6447(~ p c ), and p = 0.75. 

As might be expected, the results for p < p c and p ~ p c approach the 
values 2/3 and 0.63 respectively, further supporting that for p < p C) v = 2/3, 
while for p = p c , v ~ 0.63. The results for p > p c are inconclusive, as the 
local exponent does not yet converge at this triangle size. But, for p > p c it 
is expected that the situation is governed by one big percolation cluster. If 
this is the case, it is possible (and much less compuer power demanding) to 
"grow" percolation clusters and to find the dependence of D r on t in that 
way. The local exponents of D r for directed polymer at p — 0.75 and for 
directed percolation at the same probability are shown in figure 4. 

As can be seen in figure 4, the results of directed polymer and directed 
percolation are statistically indistinguishable, and it is quite safe to assume 
that this situation will not change for larger t. The local exponent for directed 
percolation approches a value ~ 2/3, a result which was obtained earlier by 
Balents and Kardar [8]. Thus we conclude that for p > p c , v of directed 
polymer has the same value as for p < p c , i.e. v ~ 2/3. 

It was mentioned above that the results of Figures 1,2 can be explained in 
terms of the longitudinal correlation length £y. More direct result follows from 
a picture of the optimal path as a series of zero sum segments whose mean 
length is L, connected by single bonds of value 1. (Obviously, this picture 
holds only for t >> L). According to this picture,we expect that E ~ t/L, 
and for fixed t, E x L ~ constant. Figure 5 presents the results for E/E , 
L/Lq, and (E x L)/(E x L ) in the range of probabilities 0.5 < p < 0.64, 
where E is E at p = 0.5 and L is L at p = 0.5. 

As can be seen in Figure 5, E and L form mirror reflection of each other 
over two orders of magnitude, while E x L is a slightly decreasing function 
of p. So there is almost one to one correspondence between results obtained 
from the directed polymer model (the values of E), and results obtained 
independently from the directed percolation model (the values of L). 
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In conclusion, it was shown that the bimodal distribution directed poly- 
mer model can be characterized in terms of the directed percolation model, 
and that the percolation thresehold probability p c plays a critical role in the 
directed polymer case. 
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Figure 3: Plot of the local exponent of D r as a function of t for several values 
of p. 
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Figure 4: Plot of the local exponent of D r as a function of t for directed 
polymer and directed percolation at p — 0.75. 
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Figure 5: Plot of E/E , L/L , and (E x L)/(E x L ) as a function of the 
probability p. E and L are the values at p = 0.5. 
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